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Several puzzling aspects of interplay of the experimental lattice distortion and the the magnetic 
properties of four narrow f2 9 -band perovskite oxides (YTiOa, LaTiOa, YVO3, and LaVOs) are 
clarified using results of first- principles electronic structure calculations. First, we derive param- 
eters of the effective Hubbard-type Hamiltonian for the isolated ti g bands using newly developed 
downfolding method for the kinetic-energy part and a hybrid approach, based on the combina- 
tion of the random-phase approximation and the constraint local-density approximation, for the 
screened Coulomb interaction part. Apart form the above-mentioned approximation, the procedure 
of constructing the model Hamiltonian is totally parameter-free. The results are discussed in terms 
of the Wannier functions localized around transition-metal sites. The obtained Hamiltonian was 
solved using a number of techniques, including the mean-field Hartree-Fock (HF) approximation, 
the second-order perturbation theory for the correlation energy, and a variational superexchange 
theory, which takes into account the multiplet structure of the atomic states. We argue that the 
crystal distortion has a profound effect not only on the values of the crystal-field (CF) splitting, 
but also on the behavior of transfer integrals and even the screened Coulomb interactions. Even 
though the CF splitting is not particularly large to quench the orbital degrees of freedom, the crys- 
tal distortion imposes a severe constraint on the form of the possible orbital states, which favor 
the formation of the experimentally observed magnetic structures in YTi03, YVO3, and LaVOs 
even at the level of mean-field HF approximation. It is remarkable that for all three compounds, 
the main results of all-electron calculations can be successfully reproduced in our minimal model 
derived for the isolated t2 g bands. We confirm that such an agreement is only possible when the 
nonsphericity of the Madelung potential is explicitly included into the model. Beyond the HF ap- 
proximation, the correlations effects systematically improve the agreement with the experimental 
data. Using the same type of approximations we could not reproduce the correct magnetic ground 
state of LaTiOs. However, we expect that the situation may change by systematically improving 
the level of approximations for dealing with the correlation effects. 

PACS numbers: 71.28.+d; 75.25.+z; 71.15.-m; 71.10.-w 



I. INTRODUCTION 

The transition-metal perovskite oxides ABO3 (with A— Y, La, or other trivalent rare-earth ion, and B= Ti or 
V) are regarded as some of the key materials for understanding the strong coupling among spin, orbital, and lattice 
degrees of freedom in correlated electron systems i 

According to the electronic structure calculations in the local-density approximation (LDA), all these compounds 
can be classified as u t2 g systems" , as all of them have a common transition-metal i2g-band, located near the Fermi 
level, which is well separated from the oxygen- 2p band and a hybrid transition-metal e g and either Y(4c?) or La(5c?) 
band, located correspondingly in the lower- and upper-part of the spectrum (Fig. 0. The number of electrons that 
are donated by each Ti and V site into the i2g-band is correspondingly one and two. These electrons are subjected 
to the strong intraatomic Coulomb repulsion, which is not properly treated by LDA and requires some considerable 
improvement of this approximation, which currently processes in the direction of merging LDA with various model 
approaches for the strongly-correlated systemsi2iiM£ Nevertheless, LDA continues play an important role for these 
systems as it naturally incorporates into the model analysis the effects of the lattice distortion, and does it formally 
without any adjustable parameters. Although the origin of the lattice distortion in the ti g perovskite oxides is not 
fully understood, is is definitely strong and exhibits an appreciable material-dependence, which can be seen even 
visually in Fig. El The interplay of this lattice distortion with the Coulomb correlations seems to be the key factor for 
understanding the large variation of the magnetic properties among the ti g perovskite oxides. The difference exists 
not only between Ti- and V-based compounds, but also within each group of formally isovalent materials, depending 
on whether it is composed of the Y or La atoms. The latter seems to be a clear experimental manifestation of the 
distortion effect, which is related with the difference of the ionic radii of Y and La. All together this leads to the 
famous phase diagram of the distorted ti g perovskite oxides, where each material exhibits quite a distinct magnetic 
behavior: YTiOs is a ferromagnetj 6 ! 7 ! 8 ! 9 : 10 LaTi0 3 is a three-dimensional (G-type) antiferromagnet>ii*i£ YVO3 has 
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FIG. 1: (Color online) Total and partial densities of states for YT1O3, LaTiOs, YVO3 (orthorhombic phase), and LaVOs in 
the local-density approximation. The shaded area shows the contributions of the transition-metal 3d-states. Other symbols 
show positions of the main bands. The Fermi level is at zero energy. 




FIG. 2: (Color online) Crystal structure of orthorhombic LaTiOs (left), YTi03 (middle), and YVO3 (right). La and Y sites 
are indicated by big blue spheres, Ti and V sites - by medium red spheres, and oxygen sites - by small green spheres. The 
symbols a, b, and c stand for the orthorhombic translations. The symbols 1-4 denote the transition-metal sites, which forms 
the unit cell of the distorted perovskite oxides. 



the low-temperature G-type antiferromagnetic (AFM) phase, which at around 77 K transforms into a chain-like 
(C-type) antiferromagnetic phase^iii^i^i^ and LaVC>3 is the C-type antiferromagnetii2ii£ 

On the theoretical side, the large variety of these magnetic phases has been intensively studied using model ap- 
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.roaches (Refs. [T9II20II21II2I2I2I as well as the first-principles electronic structure calculations (Refs. I2I2I27II2I2I 



I30ll3lll32ll33r) . The problem is still far from being understood, and remains to be the subject of numerous contradic- 
tions and debates. Surprisingly that at present there is no clear consensus not only between model and first-principles 
electronic structure communities, but also between researchers working in each of these groups. Presumably, the 
most striking example is LaTi03, where in order to explain the experimentally observed G-type AFM ground state, 
two different models, which practically exclude each other, have been proposed. One is the model of orbital liquid, 
which implies the degeneracy of the atomic ti g levels in the crystalline environment. 20 Another model is based on the 
theory of crystal-field (CF) splitting, which lifts the orbital degeneracy and leads to the one particular type the orbital 
ordering compatible with the G-type antiferromagnetismiSiSi The situation in the area of first-principles electronic 
structure calculations is controversial as well. Although majority of the researchers now agree that in order to describe 
properly the electronic structure of t^ g perovskite oxides, one should go beyond the conventional LDA and incorporate 
the effects of intraatomic Coulomb correlations, this merging is typically done in a semi-empirical way, as it relies on a 
certain number of adjustable parameters, postulates, and the form of the basis functions used for the implementation 
of various corrections on the top of LDA i' 2 '34 There are also certain differences regarding both the definition and the 
approximations used for the CF splitting in the electronic structure calculations, which will be considered in details 
in Sec. If 11 A1 Since the magnetic properties oiti g perovskite oxides are extremely sensitive to all such details, it is not 
surprising that there is a substantial variation in the results of first-principles calculations, which sometimes yield even 
qualitatively different conclusions about the CF splitting and the magnetic structure of the distorted t2 g perovskite 
oxides i^SiSii^SiS These discrepancies put forward a severe demand on the creation of a really parameter-free scheme 
of electronic structure calculations for the strongly-correlated systems. 

Therefore, the main motivation of the present work is twofold. 

(i) In our previous work (Ref. I34T) we have proposed a method of construction of the effective Hubbard-type model 
for the electronic states near the Fermi level on the basis of first-principles electronic structure calculations. In the 
present work we apply this strategy to the ti g states of the distorted perovskite oxides. Namely, we will derive 
the parameters of the Hubbard Hamiltonian for the ti g bands and solve this Hamiltonian using several different 
techniques, including the Hartree-Fock (HF) approximation, the perturbation theory for the correlation energy, and 
the the theory of superexchange interactions taking into account the effects of the multiplet structure of the atomic 
states. Of course, our method is based on a number of approximation, which have been introduced in Ref. |3 an d W1U 
be briefly discussed in Sec. II I II However, we would like to emphasize from the very beginning that our policy here 
is not to use any adjustable parameters apart from the approximations considered in Ref. Thus, we believe that 
it poses a severe test for the proposed method, and the obtained results should clearly demonstrate that our general 
strategy, which can be expressed by the formula first-principles electronic structure calculations — > construction of the 
model Hamiltonian — > solution of the model Hamiltonian^^ is indeed very promising. For example, at the HF level, 
using relatively simple model Hamiltonian, which is limited exclusively by the ti g bands, we will be able to reproduce 
the main results of all-electron calculations^^ Furthermore, due to the simplicity of the model Hamiltonian we can 
easily go beyond the HF approximation and include the correlation effects. 

(ii) Why do we need to convert the results of first-principles electronic structure calculations into a model? Apart 
from the purely computational reason, related with the reduced dimensionality of the Hilbert space for the solution 
of the many-electron problem,- the story of distorted ti g perovskite oxides clearly shows that the model consideration 
has yet another advantage, which is typically not sufficiently appreciated in the computational community. It is 
true that the field of first-principles electronic structure calculations is currently on the rise, and the calculation of 
the basic properties for many materials will soon become a matter of routine. However, the methods of electronic 
structure calculations are based on some approximations, the limitations of which should be clearly understood. 
Furthermore, like the experiment data, the results of first-principles electronic structure calculations will always 
require some interpretation, which would transform the world of numbers and trends into a "parallel world" of 
rationalized model categories capturing the essence of the electronic structure calculations. The understanding of 
the results of calculations in terms of these categories opens a way, on the one hand, to the material engineering of 
compounds with a desired set of properties, and, on the other way, to "engineering" of the new methods of electronic 
structure calculations in the direction of elucidation and overcoming the existing approximations. In this work we will 
illustrate how the results of first-principles calculations for the distorted perovskite oxides can be interpreted in terms 
of a limited number of model parameters, such as the crystal-field splitting, transfer integrals, and the intraatomic 
Coulomb interactions, which can be regarded as the basic operating blocks for understanding the properties of these 
materials as well as the limitation of approximations existing in the methods of electronic structure calculations. 
Particularly, we will explicitly show that the atomic-spheres-approximation (ASA), which was employed in the series 
of publications CRefs. l^lsHsils^ , is not enough as it neglects the nonsphericity of the Madelung potential. The latter 
plays an important role and in many cases predetermines the character of the magnetic ordering in the distorted ti g 
perovskite oxides. We will also show that once the parameters of Coulomb interactions are determined from the 
first principles, the commonly used mean-field HF approximation does not necessary guaranty the right answer for 
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the magnetic properties of t<i g perovskite oxides. However, we will argue that this is a normal situation, and in the 
majority of cases, a better agreement with the experimental data can be obtained by systematically including the 
correlation effects beyond the HF approximation. In this sense, our strategy is completely different from conventional 
LDA+U calculations, where the on-site Coulomb interaction U is typically treated as an adjustable parameter (e.g., 
Refs. 1271.29 31 32 33). By changing U, one can certainly get a better numerical agreement with some experimental 
data already at the HF level. However, one should clearly understand that such an empirical treatment actually 
disguises the actual role played by the correlation effects in the narrow-band compounds. 

The rest of the paper is organized as follows. In Sec. [H] we will briefly remind the main details of the crystal 
and magnetic structure of the distorted perovskite oxides. The procedure of constructing the model Hamiltonian as 
well as the results of calculations of the CF splitting, transfer integrals, and on-site Coulomb interactions for the 
isolated t2 g band will be briefly explained in Sec. IIIII Particularly, in Sec. IIII Al we will discuss highly controversial 
situation around the values of the CF splitting extracted from electronic structure calculations j££L2$L21i2& and argue 
that the main difference is caused by two factors: (i) certain arbitrariness with the choice of the Wannier functions 
for the t2 9 bands of the distorted perovskite oxides; (ii) additional approximations used for the nonspherical part of 
the crystalline potential inside atomic spheres. The methods of solution of the model Hamiltonian will be described 
Sec. II VI and the results of calculations will be presented in Sec. [3 Finally, in Sec. E] we will summarize the main 
results of our work. 



II. CRYSTAL AND MAGNETIC STRUCTURES 

The distorted perovskite oxides contain four formula units in the primitive cell. The transition-metal (B) atoms 
are located at (0,0,0) (site 1), (a/2,b/2,0) (site 2), (0,0, c/2) (site 3), and (a/2, b/2, c/2) (site 4), in terms of three 
primitive translations: a, b, and c (see Fig.|5Jl. The distortion can be either orthorhombic or monoclinic. 

The space group of the orthorhombic phase is (in the Schonflies notations or Pbnm in the Hermann-Maguin 
notations, No. 62 in the International Tables). In this case all i?-sites are equivalent and can be transformed to each 
other using symmetry operations of the D\\ group. 

The monoclinic phase has the space group C\ h {P2\/a, No. 14 in the International Tables)ii^£i In this case, there 
are two nonequivalent pairs of B-sites: (1,2) and (3,4). Each pair is allocated within the same ab-plane, so that the 
atoms can be transformed to each other using symmetry operations of the C\ h group. However, there is no symmetry 
operation, which connects the atoms belonging to different ab-planes. 

We use the experimental lattice parameters and the atomic positions reported in Ref. |(J for YTiC>3 (the data 
corresponds to the temperature T=293 K), in Ref. 12 for LaTi0 3 (T=8 K), in Ref. 14 for YV0 3 (T=65 K and 100 
K, for the orthorhombic and monoclinic phase, respectively), and in Ref. 18 for LaV03 (T=10 K). 

There are five possible magnetic structure, which can be obtained by associating with each transition-metal site 
either positive (|) or negative (|) direction of spin, without enlarging the unit cell. Therefore, each magnetic structure 
can be denoted by means of four vectors associated with the transition- metal sites (1 2 3 4). They are 

1- (TTTT); which is called the ferromagnetic (F) phase; 

2- (till), the layered (A-type) antiferromagnetic phase; 

3- (T4T4), the chainlike (C-type) antiferromagnetic phase; 

4- (Tilt); the totally antiferromagnetic (G-type) phase; 

5. (|TTT)j the spin-flip phase. In the monoclinic structure, there arc two different spin-flip phases: (ItTT) an£ l (TTIT)j 
which will be denoted as flip-I and flip-II, respectively. 

Similar classification can be used the orbital ordering. Typically, two orthogonal orbitals at the neighboring transition- 
metal sites are said to be ordered antiferromagnetically, although such a definition is not unique^ 



III. MODEL HAMILTONIAN 

Our first goal is the construction of the effective multi-orbital Hubbard model for the isolated t2 g bands: 

RR' a/3 R a(3-fS 

where c^ Q (cr q ) creates (annihilates) an electron in the Wannier orbital W§[ of the site R, and a is a joint in- 
dex, incorporating all remaining (spin and orbital) degrees of freedom. The matrix /iRR' = ||/ip^/|| parameterizes 
the kinetic energy of electrons, where the site-diagonal part (R=R/) describes the local level-splitting, caused by 
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the crystal field and (or) the spin-orbit interaction, and the off-diagonal part (R^R') stands for the transfer inte- 
grals. U af37 s = JdrJ dr'W^ t (r)^(r)w scr (r-r')W^ t (r')M4(r') = (W^W^\v scr \W^W^} are the matrix elements of 
screened Coulomb interaction v SCI (r— r'), which are supposed to be diagonal with respect to the site indices. 

The parameters of the Hubbard Hamiltonian Q can be derived "from the first principles" , starting from the 
electronic structure in LDA. This procedure has been already discussed in details in Ref. l34l Here we only remind 
the main ideas and present the results for the distorted perovskite compounds. 

All calculations have been performed using linear muffin-tin-orbital (LMTO) method in the atomic-spheres- 
approximation^ We have also considered the additional correction to the crystal-field splitting, coming from the 
nonsphericity of electron-ion interactions, beyond conventional ASAi^i 

A. Kinetic-Energy Part, Controversy about the Crystal-Field Splitting 

The kinetic-energy part of the Hubbard Hamiltonian was constructed using the downfolding methodi2Si24 It yields 
certain set of parameters {ft^,}. The Wannier functions {Wr} for the the isolated t% g bands can be formally 
reconstructed from {/irjj/} using the definition h^^j 

=(W^}H huA \W^), where H hDA is the LDA Hamiltonian in 

ASAM. 

The (characteristic) example of such Wannier functions constructed for LaTi03 is shown in Fig. [3] and their 
extension in the real space is illustrated in Fig. 21 The functions are well localized: about 80-85% of their total weight 
is concentrated at the central Ti site, 5-9 % belong to neighboring oxygen sites, and about 10 % is distributed over 
La, Ti, and O sites located in next coordination spheres. Another measure of localization of the Wannier functions is 
the expectation value of the square of the position operator: (r 2 ) Q = (H^|(r— R) 2 |W^), — which yields (r 2 ) a = 2.68, 
2.36, and 2.37 A 2 for a= 1, 2, and 3, respectively. The Wannier functions for LaTiC>3 are less localized in comparison 
with the more distorted YTiC>3, where (r 2 ) is of the order of 1.90-2.28 A 2 . 34 - However, this is to be expected. 

The parameters {h^ R ,} include all kinds of hybridization (or covalent mixing) effects between transition-metal ti g 

and other atomic orbitals. However, there are other effects, which are not yet included in {h^^,}. They come from 
the nonsphericity (n-s) of the Madelung potential for the electron-ion interactions, and contribute to the CF splitting. 
The proper correction to {/i^,} can be easily calculated in the basis of Wannier functions {W 7 ^;}: 

A"-< R = £ (^I ir^rJ ^' (2) 

R'^R 1 1 

where Z^, is the total charge associated with the site R' (namely, the nuclear charge minus the screening electronic 
charge encircled by the atomic sphere), and r is the position of the electron in the sphere R. 

The main idea behind this treatment of the CF splitting is based on certain hierarchy of interactions in solids. 
It implies that the strongest interaction, which leads to the energetic separation of the t^g band from other bands 
(Fig. UJ, is the covalent mixing. For example, in many transition- metal oxides this interaction is mainly responsible 
for the famous splitting between the transition-metal t^g and e g states^ The nonsphericity of the Madelung potential 
is considerably weaker than this splitting. However, it can be comparable with the effects of covalent mixing in the 
narrow t^g band. Therefore, the basic idea is to treat this nonsphericity as a pseudo-perturbation^ and calculate 
the matrix elements of the Madelung potential in the basis of Wannier functions constructed for spherically averaged 
ASA potential. 

The same strategy can be used for the spin-orbit (s-o) interaction, which yields the following correction to the 
kinetic part of the model Hamiltonian: 

A s -°^r = (^^(W x p) • a\WS> 

(V being the LDA potential and <x being the vector of Pauli matrices) . 

One of the most controversial issues, which is actively discussed in the literature, is the values and the directions 
of the CF splitting in the distorted ti g perovskite oxides. Therefore, we would like to discuss this problem more 
in details. Basically, there are two sources of discrepancies, which largely affect the conclusions about the orbital 
ordering and the magnetic ground state: 

1. (The origin of the CF splitting: the nonsphericity of the Madelung potential versus the hybridization effects) 
The importance of nonsphericity of the Madelung potential has been emphasized by several authors. The original 
idea is due to Mochizuki and Imada, who considered the i2g-level splitting in the Ti-compounds associated with 
the displacements of Y and La sitesi^ 3 . It was paraphrased by Cwik et al.^ who suggested the main effect comes 
from the deformation of the TiC>6 octahedra. A more general picture has been considered by Schmitz et al.fil who 
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FIG. 3: (Color online) Contour plot of the Wannier functions in the case of LaTiOs in three orthorhombic planes: ab (top), ac 
(middle), and be (bottom). The solid and dashed line correspond to the positive and negative values of the Wannier functions. 
The projections of different atoms on the planes are denoted by the following symbols: + (La), o (Ti), and x (O). Around each 
site, the Wannier function increases/decreses with the step 0.04 from the values indicated on the graph. 



summed up all contributions in the Madelung potential. A weak point of all these approaches is the approximate 
treatment of the hybridization effects, which largely relied on the model parameters. Moreover, the result depends 
crucially on the value of the dielectric constant, which is treated as an adjustable parameter. On the other hand, the 
parameters of the model Hamiltonians extracted from the first-principles electronic structure calculations using either 
downfolding (Refs. |2^, [3(], and l3lh or Wannier function (Ref. l32|) methods automatically include all effects of the 
covalent mixing. In this sense, these are more rigorous techniques. However, all these calculations were supplemented 
with the additional atomic-spheres-approximation and neglected the nonspherical part of the Madelung potential. 
The Madelung term has been also neglected in our previous work (Ref. l28l). As we shall see in Sees. IV Bl and IV Dl it 
will definitely revise several statements of Ref. |3 However, the final conclusion about the magnetic ground state of 
YTiC>3 and LaTiC>3 remains valid. 

2. (The nonuniqueness of the Wannier functions) Different calculations yield different parameters of the CF splitting. 
For example, for LaTiC>3 different authors reported the following parameters of the CF splitting (between lowest and 
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FIG. 4: Spacial extension of the Wannier functions in the case of LaTiOa: the electronic charge accumulated around the 
central Ti site after adding every new sphere of neighboring atomic sites. 

highest ^-levels): 93 meV (Ref-H), 200 meV (Ref.HJ), and 270 meV (Ref.©. Th ere is a particularly bad custom 
to criticize Ref. I^sl ^^ 2 ^ which reports the smallest value, even with certain hints at the accuracy of calculations. 32 
It is also premature to think that the small CF splitting will inevitably lead to the orbital liquid scenario ^fl* 3 ^ because 
the transfer interactions are also strongly affected by the lattice distortion, which makes a big difference from the 
idealized cubic perovskitesi^i 

First, we would like to consider the second part of problem and argue that different values of the CF splitting are 
most likely related with the different choice of the Wannier functions, which by no means is unique. This is not a 
problem of accuracy of calculations. 

In the downfolding method employed in Refs.l23.l30l andl3ll all basis states were divided in two groups: the "i2g" 
part {x*}) and the "rest" of the basis functions {Xr}- The effective Hamiltonian is constructed by eliminating the 
"rest" parti^i A similar idea (although formulated in slightly different way) is employed in projections scheme for the 
construction of the Wannier functions j 3 ^ where {x*} play a role of trial orbitals. The basic difficulty here is that, 
in the distorted perovskites, the set of atomic "^29" orbitals cannot be defined in an unique fashion: since the local 
symmetry is not cubic, the abbreviations like "^29" an d "eg" will always reflect some bad quantum numbers for the 
state, which are mixed by the crystal field and/or the transfer interactions. In numerical calculations, the set of "tig" 
orbitals is always specified in some local coordinate frame, and the choice of this frame appeared to be different in 
different calculations. For example, Pavarini et al. (Refs. l30l and l3lf) and Streltsov et al. fRef. f32j) selected their local 
coordinate frames from some geometrical considerations. A completely different strategy was pursued by the present 
author in Ref.l28t where the atomic "£29" orbitals were determined from the diagonalization of the local density matrix 
constructed from the t2 g bands and projected onto all five 3d-orbitals of the transition-metal sites. 

Then, we are ready to argue that different choice of the local coordinate frame naturally explains the difference 
in the parameters of the CF splitting reported by different authors. For these purposes we consider two different 
setups in the downfolding scheme, and consider LaTiOa as an example. The first scheme is absolutely identical to 
that proposed in Ref.l28l where the atomic "i2g" orbitals have been defined as three most populated orbitals obtained 
from the diagonalization of the density matrix for the ti g bands. In the second scheme, we first construct a more 
general 40x40 tight-binding Hamiltonian, comprising the Ti(3d) and La(5d) states, and reproducing the behavior of 
40 overlapping Ti(3d)-La(5d) bands. Other orbitals have been eliminated using the downfolding method. Then, wc 
diagonalize the site-diagonal part of this Hamiltonian and define three lowest eigenstates at the Ti-site as the atomic 
"^29" orbitals. After that we eliminate the rest of the Ti(3ci) and La(5ci) states using the downfolding method and 
obtain the minimal 12x12 Hamiltonian for the ti g bands. 

Both downfolding schemes are nearly perfect and well reproduce the behavior of the t2 g bands in the recipro- 
cal space (Fig. |SJ). However, they yield very different parameters after the Fourier transformation to the real 
space. For example, the splitting of atomic ti g levels (in meV) obtained in the schemes I and II is (—49,4,44) 
and (—320,123,197), respectively. Moreover, the eigenvectors corresponding to the lowest "£29" levels appear to 
be also different. In the orthorhombic coordinate frame, specified by the vectors a, b, and c, the eigenvec- 
tors have the following form (referred to the site 1): |*i)=0.32|xy)-0.73|yz)-0.10|z 2 )-0.18|zx)+0.57|x 2 -y 2 ) and 
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FIG. 5: (Color online) Left panel: LDA energy bands for LaTiOs obtained in LMTO calculations and after the tight-binding 
(TB) parametrization using two different schemes of downfolding. In the scheme I, the local coordinate frame has been obtained 
from the diagonalization of the local density matrix. In the scheme II, the local coordinate frame has been obtained from the 
diagonalization of the site-diagonal part of a more general 40x40 Ti(3d)-La(5d) tight-binding Hamiltonian. Notations of the 
high-symmetry points of the Brillouin zone are taken from Ref . I42L Right panel: distance-dependence of averaged parameters of 

the kinetic energy hnw (d) = (y2 a ^ ^rr'^r/r) » where d i s the distance between transition- metal sites R and R'. The data 

for d—Q correspond to the crystal-field splitting of the covalent type, and d~4 A is the distance between nearest transition-metal 
sites. Two schemes generate very similar electronic structure in the reciprocal space, which is well consistent with results of 
original LMTO calculations. However, the transfer integrals and the crystal-field splitting are different. Scheme II yields larger 
crystal-field splitting, but less localized transfer integrals. 



|*n)=0.31|xy)-0.20|yz)+0.15|z 2 )+0.55|zx)+0.73|x 2 -y 2 ), for the scheme I and II, respectively. Thus, the small value 
of the CF splitting reported in Ref.|28|is related with the particular choice of the Wannier functions (or the parameters 
of the downfolding scheme). Had we changed our definition of the local coordinate frame, our conclusion would have 
been also different, and we could easily obtain the CF splitting of the order of 500 meV (and even larger). 

Then, it is of course right to ask which scheme is better? In principle, the physics should not depend on the choice 
of the Wannier functions, and as long as we are dealing only with the kinetic-energy part of the model Hamiltonian, 
both schemes are totally equivalent as they equally well reproduce the behavior of the ti g bands. However, what we 
want to do next is to combine this kinetic- energy part with the Coulomb interactions, and to use only the site-diagonal 
part of these interactions. This is of course an approximation, and in order to justify it one should guarantee that the 
Wannier functions, which are used as the basis for the matrix elements of the Coulomb interactions, were sufficiently 
well localized in the real space, so that all inter-site matrix elements could be neglected. 

The degree of localization of the Wannier orbitals is related with the spread of transfer integrals. Loosely speaking, 
in order to contribute to the transfer integral between N-th neighbors, the Wannier function should have a finite 
weight at this neighbor. The distance-dependence of transfer integrals calculated in two different scheme is shown in 
Fig. [3] One can clearly see that the transfer integrals obtained in the scheme I are indeed well localized and basically 
restricted by the nearest neighbors. The transfer integrals obtained in the scheme II are less localized and spread 
far beyond the nearest neighbors. Therefore, the Wannier functions, corresponding to the scheme I, should be more 
localized. This is not surprising, because the scheme I guarantee that the density matrix (or the integrated density of 
states) at the transition-metal site is already well described by the central parts (or "heads") of the Wannier functions, 
given by the atomic orbitals {xt}~ Therefore, the tails of the Wannier functions, coming from the neighboring sites, 
should be small and cancel each other. In the scheme II, the local density matrix is composed by both "heads" of the 
Wannier functions as well as their tails coming from the neighboring sites. Intuitively, this means that the tail part 
of the Wannier functions is larger for the scheme II and these functions are less localized. 

Thus, we believe that the scheme I is more suitable for the purposes of our work and we apply it to all perovskite 
compounds. The transfer integrals, obtained in such a way, are indeed well localized and restricted by the nearest 
neighbors (Fig. 0). 

The next important contribution to the CF splitting comes from the nonsphericity of the Madelung potential. This 
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FIG. 6: (Color online) Behavior of averaged parameters of the kinetic energy for various perovskite compounds. In the 
orthorhombic (o) structure, all four sublattices of the transition-metal sites are equivalent as they can be transformed to each 
other by the symmetry operations of the Di,h group. In the monoclinic (m) structure, the transfer interactions corresponding 
to two different sublattices are shown by closed and open symbols. For other notations, see Fig. [21 
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FIG. 7: Convergence of the crystal-field splitting as the function of cut-off radius in Eq. |(2J. 



effect is comparable with the CF splitting of the covalent type. We have also found that there are several different 
contributions to the CF splitting, which tend to cancel each other. For example, the CF splitting of the covalent 
type in YTi03 and LaTi03 is largely compensated by the nonsphericity of the Madelung potential coming from the 
region encircling the the neighboring oxygen sites and corresponding to the cut-off radius <i~2 A in Fig.[7[ 43 The next 
important contribution comes from the Y/La and Ti sites, located in the next coordination spheres (d~4 A). In some 
compounds, even longer-range interactions spreading up to d^lO A can have a relatively large weight in the sum J3J). 

The final scheme of the "t2 fl " -level splitting, which takes into account all these effects, is summarized in Fig. [HI The 
splitting is not particularly large. Nevertheless, as we shall see below, it appears to be sufficient in order to explain the 
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FIG. 8: t2 S -level splitting in various compounds. The notations "site 1" and "site 3" correspond to two nonequivalent 
transition-metal sites in the monoclinic structure (shown in Fig. HJl. 



experimentally observed orbital ordering and the magnetic ground state in all considered compounds, except LaTiO^. 
The Madelung potential changes not only the magnitude of the crystal-field splitting, but also the type of the atomic 
"t2 ff "-orbitals, which are split off by the distortion (Table^J. The type of these orbitals is extremely important in the 



TABLE I: Lowest (for YTiC>3 and LaTiOa) and highest (for YVO3 and LaVOa) atomic orbital obtained from the diagonalization 
of the site-diagonal part of the model Hamiltonian in the downfolding method (denoted as the "covalent part") and and after 
including the nonsphericity of the Madelung potential (denoted as "total"). The symbols 'o' and 'm' stand for orthorhombic 
and monoclinic phases, respectively. The order of the basis orbitals is xy, yz, z 2 , zx, and x 2 — y 2 , in the orthorhombic coordinate 
frame. The positions of the transition-metal sites are shown in Fig. EI 



compound 


phase 


site 




covalent part 






total 






YTiOs 





1 


( 0.32, 


0.78, 0.36,-0.33,- 


-0.22) 


(-0.13, 


0.45, 0.38,- 


-0.61, 


0.50) 


LaTi0 3 





1 


(-0.32, 


0.73, 0.10, 0.18,- 


-0.57) 


(-0.06, 


0.85, 0.15, 


0.34, 


0.37) 


YVO3 





1 


( 0.04, 


0.84, 0.19,-0.22, 


0.47) 


( 0.29, 


0.90, 0.25,- 


-0.22, 


-0.06) 


YVO3 


111 


1 


(-0.04, 


0.47, 0.19,-0.61, 


0.60) 


( 0.12, 


0.54, 0.18,- 


-0.74, 


0.33) 






3 


(-0.01, 


0.04,-0.09, 0.79, 


0.60) 


( 0.28, 


0.24,-0.28, 


0.89, 


-0.05) 


LaV0 3 


111 


1 


( 0.04,- 


-0.31, 0.11, 0.80,- 


-0.49) 


( 0.16, 


0.09, 0.01, 


0.98, 


-0.04) 






3 


(-0.01, 


0.37, 0.10, 0.75, 


0.54) 


(-0.09, 


0.55, 0.14, 


0.78, 


0.24) 



analysis of the magnetic properties. 

Finally, both magnitude and form of the CF splitting can be different for the nonequivalent transition-metal sites in 
the monoclinic structure. Therefore, one can generally expect rather different magnetic behavior in two nonequivalent 
ab-planes of the monoclinic phase. 



B. Effective Coulomb Interactions 

The effective Coulomb interaction in the ti g band is defined as the energy cost for moving an electron between two 
Wannier orbitals, W£ and Wi,, , which have been initially populated by n Rct and riR/,3 electrons: 

U aa pp = E [n Ra + 1, n R >j3 — 1] — E [n Ra , nwp] ■ 

More generally, one can consider the electron transfer from any linear combination of the Wannier orbitals at the site 
R' to any linear combination at the site R. This define the full matrix of screened Coulomb interactions f/=||C/a/3 7 5 li- 
lt can be calculated under certain approximations, which have been discussed in details in Rcf. 34. The method 
consists of two parts. 

1. First, we perform the standard constraint-LDA (c-LDA) calculations and artificially switch off all matrix elements 
of hybridization involving the atomic 3d-states4s This part takes into account the screening of Coulomb interactions 
caused by the relaxation of the 3d-atomic basis functions and the redistribution of the rest of the charge density. 
Typical values of on-site Coulomb interactions (it) obtained in this approach for the Ti- and V-based perovskites vary 
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between 8.5 and 9.3 eV. Using a similar approach, one can calculate the intraatomic exchange coupling constant (J), 
which is about 0.9 eV for all considered compounds. Finally, from the obtained parameters u and j, one can restore 
the full 5x5x5x5 matrix u=||u a| a 7 ,5|| of screened Coulomb interactions between atomic 3<i-orbitals with the same 
spin, as it is typically done in the LDA+U method^ 

2. Then, we switch on the hybridization and evaluate the screening caused by the change of this hybridization between 
the atomic 3d-orbitals and the rest of the basis states in the random-phase approximation (RPA) : 

U= l-uP(O) u. (3) 

This scheme implies that different channels of screening can be included consecutively. Namely, the u-matrix derived 
from c-LDA is used as the bare Coulomb interaction in the Dyson equation and the 5x5x5x5 polarization matrix 
P=||P Q/ 3 7 j|| describes solely the effects of hybridization of the 3d-states with 0(2p), and either Y(4d) or La(5d) states, 
which lead to the formation of the distinct oxygen- 2p, transition-metal t2 g , and a hybrid e g band in Fig. ^ The 
matrix elements P are given by 

P a , j5 (u) — — + £nk + lS{nnk _ nn , k ,) . (4) 

nk n'k' 

where {e„k} an d {n„k} are LDA eigenvalues and occupation numbers for the band n and momentum k in the first 
Brillouin zone (the spin index is already included in the definition of n) and d- yn ] i =(x 7 \ip n \ i } is the projection of LDA 
eigenstate i/Vik onto the atomic 3d orbital Xy 

Hence, we obtain a 5x5x5x5 matrix of screened Coulomb interactions in the basis of all five 3d orbitals. This 
matrix is then transformed into the local coordinate frame spanned by three "^23" orbitals with the same spin. The 
obtained 3x3x3x3 matrix is expanded in the spin subspace using the orthogonality condition between Wannier 
orbitals with different spins. This yields a 6x6x6x6 matrix U, which is used in the actual calculations. 

Only for explanatory purposes, we fit U in terms of two Kanamori parameters: the intra-orbital Coulomb interaction 
U and the exchange interaction Jm^ The results of such fitting are shown in Table [H] There is the clear dependence 
of the parameter U on the local environment in solid, which is captured by the RPA calculations. — Generally, the 
value of U is larger for more distorted Y-compounds. There is also a clear correlation between the value of U and 
the magnitude of local distortion around two nonequivalent transition-metal sites in the monoclinic phases: the sites 
experiencing larger distortion (according to the magnitude of the CF splitting in Fig. (BJ have larger U. On the other 
hand, J is practically insensitive to the local environment in solids. 



TABLE II: Results of fitting of the effective Coulomb interactions in the ti g band obtained in the hybrid c-LDA+RPA approach 



in terms of two Kanamori parameters: the intra-orbital Coulomb interaction U and the exchange interaction J (in eV) . 45 The 
symbols 'o' and 'm' stand for the orthorhombic and monoclinic phases, respectively. The positions of the transition-metal sites 
are shown in Fig. [5] 

compound phase site U J 

YTiOs o 1 3.45 0.62 

LaTi0 3 o 1 3.20 0.61 

YV0 3 o 1 3.27 0.63 

YVO3 m 1 3.19 0.63 

3 3.26 0.63 

LaVOa m 1 3.11 0.62 

3 3.12 0.62 



IV. SOLUTION OF MODEL HAMILTONIAN 



In this section we briefly discuss the methods of solution of the model Hamiltonian (JIJ. We start with the simplest 
HF approach, which totally neglect the correlation effects. Then, we consider two simple corrections to the HP 
approximation, which allow to include some of these effects. One is the second-order perturbation theory for the total 
energy. It shares common problems of the regular (nondegenerate) perturbation theory and allows to calculate easily 
the correction to the total energy, starting from the single-Slater-determinant HF approximation for the wavefunctions. 
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Therefore, we expect this method to work well for the systems where the orbital degeneracy is already lifted and the 
ground state is described reasonably well by a single Slater determinant, so that other corrections can be treated as 
a perturbation. The second scheme is the variational superexchange theory for the d 1 perovskites, which takes into 
account the multiplet structure of the excited atomic states. It allows to study the effect of electron correlations on the 
orbital ordering. However, it is limited by typical approximations made in the theory of superexchange interactions, 
which treat all transfer integrals as a perturbation. 

All calculations have been performed in the basis of Wannier functions {W 7 ^;}, which have a finite weight at the 
central transition-metal sites as well as the oxygen and other atomic sites located in the nearest neighborhood to the 
transition-metal atoms. In order to calculate the local quantities, associated with the transition-metal atoms, such 
as the spin or orbital magnetic moments as well as the distributions of the charge densities, the Wannier functions 
have been expanded over the standard LMTO basis, and the aforementioned quantities have been obtained by the 
integration over the atomic spheres of the transition-metal sites. 



A. Hartree-Fock Approximation 



The Hartree-Fock approximation provides the simplest solution of the many-body problem described by the model 
Hamiltonian ft}. In this case, the trial wavefunction for the many-electron ground state is searched in the form 
of a single Slater determinant |<S{<£> n k}), which is constructed from the one-electron orbitals {ip n k}- The latter are 
subjected to the variational principle and requested to minimize the total energy 

i?HF = min (S{tp n k}\Tl\S{(pnk}) 

for a given number of particles AT, yielding the set of well-known HF equations: 

(hk + V) Ifnk) = Enkl^nk), (5) 

where /ik=||/i k ' ,3 || * s * ne kinetic part of the model Hamiltonian in the reciprocal space: h%f=j? X)r/ ^RR/ e_lk '^ R_R ^ 
(N being the number of sites), and 1^=11 V^gH is the HF potential: 

V a f3 — {UaP^S — UaS-fp) U^S- (6) 

Eq. J^J is solved self-consistently together with the equation 

occ 
nk 



for the density matrix n=||n a| g|| in the basis of Wannier orbitals. 
After the self-consistency, the total energy can be calculated as 

occ _^ 
EhF — ^ £ nk — - ^ Vp a n a 0. 
nk af3 

Using {e ra k} and {<p n k}, one can calculate the Green function in the HF approximation: 

gRR'M = E IgSjE^j e *-(R-R'). 

10 — £ n k + 10 

nk 

The latter is widely used for the analysis of interatomic magnetic interactions^ 

Jrr' = ^ Im / £F &Tl -i {Srr'M a ^r'rM a ^} . (7) 

where £^, = 5Trs{(i±<7 z )C? RR /} is the projection of the Green function onto the majority (|) and minority (J.) spin 
states, AV—Trs{<7 z V} is the magnetic (spin) part of the HF potential, Trs (Tr^) denotes the trace over the spin 
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(orbital) indices, 1 and a z are correspondingly unity and Pauli matrices of the dimension 6, and Ef is the Fermi 
energy. 

The interatomic magnetic interactions {Jrr/} characterize the spin stiffness of the magnetic phase. Therefore, 
they can be directly compared with the experimental spin-wave spectra derived from the inelastic neutron scattering 
measurement. 

According to Eq. [7| Jrr/>0 (<0) means that for a given magnetic structure, the spin arrangement in the bond 
(RR') corresponds to the local minimum (maximum) of the total energy. However, in the following we will use the 
universal notations, according to which Jrr'>0 and <0 will stand the ferromagnetic and antiferromagnetic coupling, 
respectively. This corresponds to the local mapping onto the Heisenberg model of the form 

-E^eis = — 2^ ^ RR ' eR ' eR '' (8) 

(RR') 

where eR is the direction of the spin magnetic moment at the site R. The "local mapping" means that, strictly 
speaking, Eq. [HI is justified only for the infinitesimal rotations of the spin magnetic moment near an equilibrium. 

The magnetic interactions are extremely sensitive to the orbital ordering. Therefore, they can be regarded as the 
local prob of the orbital ordering in each magnetic state. 



B. Second Order Perturbation Theory for Correlation Energy 



The simplest way to go beyond the HF approximation is to include the correlation interactions in the second 
order of perturbation theory for the total energy^ 7 " The correlation interaction (or a fluctuation) is defined as the 
difference between true many-body Hamiltonian (I), and its one-electron counterpart, obtained at the level of HF 
approximation: 



R \ a/3fS 



Hc = 1^ | 9 2^ ^W^cS^^aS " F «/3 a Ra 2 R/3 | ■ ( 9 ) 

a/3 



By treating He as a perturbation, the correlation energy can be easily estimated as: 47 

(G\H c \S)(S\Hc\G] 
Ehf(S) — £hf(G) 



^ (G\Hc\S)(S\H c \G) 



s 

where \G) and \S) are the Slater determinants corresponding to the low-energy ground state in the HF approximation, 
and the excited state, respectively. Due to the variational properties of the Hartree-Fock approach, the only processes 
which contribute to Eq correspond to the two-particle excitations, for which each \S) is obtained from \G) by replacing 
two one-electron orbitals, say ifmki and <y9n 2 k 2 i from the occupied part of the spectrum by two unoccupied orbitals, 
say </?ri 3 k 3 and <v9n 4 k 4 - Hence, using the notations of Sec. IIIII the matrix elements take the form: 

(S\H.c\G) = (</2ri3k 3 ( / ? 7i4k 4 |Wscr| ( /?7iiki ( /9n2k2) — (V?n3k 3 <v9n 4 k 4 l^scr | <fin 2 k2 Vniki )• (11) 

These matrix elements satisfy the following condition: (5|7ic|G)~^ J2n e' i ( k3+k4 ~ kl ~ k2 )' R , provided that the screened 
Coulomb interactions are diagonal with respect to the site indices. In the following we will retain only the R=0 part 
in this sum. This corresponds to the single-site approximation for the correlation interactions, which is known to be 
good for three-dimensional systems and becomes exact in the limit of infinite spacial dimensions^ 

Finally, we employ a common approximation of noninteracting quasiparticles and replace the denominator of 
Eq. I |10|l by the linear combination of HF eigenvalues: E-hf(S) — £ , HF(G)«£ ri3 k 3 +£n 4 k 4 — £mki — £n 2 k 2 ~ 

The form of Eq. I|lt)fl implies that the HF ground state is nondegenerate, and the correlation effects can be system- 
atically included by considering the regular perturbation theory expansion. It does not apply to the cubic systems, 
where the ground state is infinitely degenerate (with respect to different orbital configurations), and where a suit- 
able approach for the correlation energy should be based on the degenerate perturbation theory^ Thus, the use of 
Eq. nitty implies that the orbital degeneracy is already lifted by the crystal distortion. As we shall see below, this 
approximation can be justified for a number of systems. 
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FIG. 9: The multiplet structure of the excited atomic configuration t|g m YTiOa and LaTiO 



C. Effects of Multiplet Structure in the Theory of Superexchange Interactions 

Another method, which allows to treat some correlation effects beyond the mean-field HF approximation is based 
on the generalization of the theory of superexchange interactions in order to describe the correct multiplet structure of 
the atomic states. Similar idea has been discussed in the context of colossal magnetoresistive perovskite manganites4^ 
The formulation is extremely simple for the d 1 compounds, like YTiC>3 and LaTiC>3. 

The superexchange interaction in the bond (RR') is basically the gain of the kinetic energy, which is acquired by an 
electron occupying the atomic orbital 0r of site R in the process of virtual hoppings into the subspace of unoccupied 
orbitals of the (neighboring) site R', and vice versai£2*£i In the atomic limit for the d 1 compounds, there is only one 
tig electron at each Ti site. This is essentially an one-electron problem, where the form of the atomic orbital </>r is 
determined by the site-diagonal part of kinetic-energy ||/IrrII> incorporating the effects of the spin-orbit interaction 
and the CF splitting. Therefore, in the pure atomic limit, the ground-state wavefunction for each bond (RR') can be 
described by a single Slater determinant: 



|Grr') = ^={</>r(1)0R<(2) 



<M1)0r(2)}. 



Then, 0r and 0r< are expanded over the Wannier orbitals associated with the site R and R' and the transfer integrals 
connecting different Wannier orbitals are treated as a perturbation. The excited states at the sites R and R', which 
appear in the process of the virtual hoppings are the two-electron states and subjected to the multiplet splitting. 
This is exactly the point where the electron correlations, beyond the HF approximation, enter the problem. In order 
to incorporate these effects we note that from m=6 Wannier spin-orbitals {W 7 ^} at each Ti site, one can construct 
■^m{m— 1)=15 antisymmetric two-electron Slater's determinants {IS 1 )} (5=1,. . . ,15), which can be used as the basis 
for the screened Coulomb interactions in the excited state: Uss'=(S\v S ci\S'). The diagonalization of this 15x15 
matrix yield the complete set of eigenvalues {£rm} and eigenfunctions {|RM)} of the two-electron states at the site 
R (M=l,. . . ,15). An example of such a multiplet structure for YTiC>3 and LaTiC>3 is shown in Fig. EI According to 
the first Hund rule, the lowest-energy configuration corresponds to the spin-triplet state 3 Ti g . The degeneracy of the 
3 Ti g , 1 E g , and 3 T2 g levels is lifted by the orthorhombic distortion, which affects the matrix elements of the effective 
U via the RPA channel of screening. The splitting is larger for the more distorted YTiC>3. 

In order to calculate the energy gain caused by the superexchange interactions, the eigenfunctions {|RM)} 
shall be projected onto the physical subspace of two-electron states which can be created by transferring an elec- 
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tron from the neighboring sites. The corresponding projector operators have the form: VR=J2 a \Pr)(Pr\> where 
I^r) = ^{^ r (^)^ / r(2) — W / r(1)0R.(2)} is the Slater determinant constructed from the occupied orbital </>r and one 
of the basis Wannier orbitals W£. In the other words, the projection Vr imposes an additional constraint, which 



guarantees that one of the orbitals in the two-electron state is 
hoppings in the bond (RR') is given by 



i>R. Then, the energy gain caused by the virtual 



A-Ei 



RR' 



G 



RR' 



J RR' 



E 

\ M 



Vr,\-R'M)(R'M\P w 



R'M 



hwR + (i «-> j; 



G 



RR' 



(12) 



The total energy of the system in the superexchange approximation is obtained after summation over all bonds, which 
should be combined with the site-diagonal elements, incorporating the effects of the CF splitting and the relativistic 
spin-orbit interaction: 



E< 

R 



»R|«RR 



E 

(RR') 



AE- 



RR' 



Finally, the set of occupied orbitals {4>r} is obtained by minimizing -Ese> e -g-> using the steepest descent method. 

For a given orbital ordering, the multiplet effects are expected to be more important in the case of the AFM spin 
ordering, where an electron comes to the neighboring site with the opposite direction of spin. In this case, the excited 
configuration is subjected to the multiplet splitting into the spin-singlet and spin-triplet states, which additionally 
stabilizes the AFM spin stated 



V. RESULTS AND DISCUSSIONS 



In this section we present results of solution the model Hamiltonian JQ) for the distorted perovskite compounds. 
We start with Y-based perovskites, where the orbital ordering is largely controlled by the CF splitting coming from 
the experimental lattice distortion. We will show that this distortion imposes a severe constraint on the magnetic 
properties and predetermines the type of the magnetic ground state. Then, we turn to the La-based perovskites, 
where the situation is less clear: while the magnetic structure of LaVC-3 can be still understood on the basis of its 
experimental lattice structure, LaTi03 poses many open and unresolved questions. 

First, we will consider the impact of the crystal structure without the spin-orbit interaction. The effects of the 
spin-orbit interaction will be discussed separately in Sec. IV El 



A. YV0 3 



YVO3 exhibits two structural phase transitions^*^ The first one is the second-order transition from orthorhombic 
(D^t) to monoclinic (apparently D^) phase, which takes places at 200 K and is believed to coincide with the onset 
of the orbital ordering. The second one is the first-order transition at 77 K from monoclinic to another orthorhombic 
phase. The magnetic transition temperature is 116 K, which lies in the monoclinic region and does not coincide with 
any structural phase transition. On the other hand, the change of the crystal structure at 77 K coincides with the 
magnetic phase transition. The magnetic structure in the interval 77 K < T < 116 K is C-type AFM, while below 77 
K it becomes G-type AFM. 



1. Low-temperature orthorhombic phase (T < 77 K) 



YVO3 has the largest CF splitting amongst perovskite compounds, which crystallize in the orthorhombic phase 
(Fig. [SJ. It lowers the energies of two ti g levels - just enough to accommodate two d electrons. The highest level is 
separated from the middle one by a 111 meV gap. Therefore, the orbital ordering in this d 2 compound is expected to 
be quenched (at least partially) by the crystal distortion. 

This is clearly seen in our Hartree-Fock calculations. The orbital ordering depends on the magnetic state. However, 
this dependence is weak and can be hardly seen on the plot (Fig. HOfl . The form of the orbital ordering, which can be 
schematically viewed as an alternation of the (xy,yz) and (xy,zx) orbitals in the cubic coordinate frame associated with 
the V06 octahedra, is in an excellent agreement with the theoretical prediction of Sawada and Terakura based on the 
semi-empirical LDA+C7 method^ which was later on confirmed by the synchrotron x-ray-diffraction measurement^ 
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FIG. 10: (Color online) Distribution of the charge density around V-sites in various magnetic phases of orthorhombically 
distorted YVO3 (T<77 K), as obtained in Hartree-Fock calculations. Different magnetic sublattices are shown by different 
colors. 



The first important question, which we would like to address here is where does this orbital ordering come from? 
In Fig. ^2 we show results of calculations obtained using three different settings for the site-diagonal part of kinetic- 
energy part of the model Hamiltonian: (i) h^ R =0 (no CF splitting); (ii) the parameters {^jm} extracted from the 
downfolding method, which takes into account only the covalent type of the CF splitting; and (iii) the parameters 
{ft-RR.} obtained in the downfolding method and corrected for the nonsphericity of the Madelung potential J5J). One 
can clearly see that in order to reproduce the correct orbital ordering, all contributions to the CF splitting appear 
to be important. Had we neglected some of these contributions, not only the orbital ordering but also the magnetic 
ground state would have been different. For example, without the Madelung term, the magnetic ground state is 
expected to be of the A-type, in clear disagreement with the experimental data. 

As it has been already discussed by other authors jiiLSLSli the C-type orbital ordering shown in Fig. 1 101 favors the 
G-type AFM spin ordering, which emerges as the ground state already at the level of HF calculations (Table ITTTf) . 
The order of the magnetic states, corresponding to the increase of the total energy, is G^C^flip— >A— >F, which is 
well consistent with results of all-electron LDA+J7 calculationsi2L22i£i 

However, the orbital ordering is not fully quenched by the crystal distortion and to certain extent can adjust the 
change of the magnetic state through the Kugel-Khomskii mechanismi^l This is seen particularly well in the behavior 
of interatomic magnetic interactions, which reveal an appreciable dependence on the magnetic state. For example, by 
going from the G-state to the F-state, the in-plane interaction J12 (J34) changes by nearly 70%, and the inter-plane 
interaction J13 (J24) changes by 25%. 

In agreement with the experimental finding^) the magnetic interactions in the G-type AFM ground state are 
nearly isotropic. However, this isotropic behavior can be easily destroyed by a small change of the orbital ordering, 
which is realized for example in other magnetic states. 



17 




FIG. 11: (Color online) Orbital ordering in the G-type antiferromagnetic phase of orthorhombically distorted YVO3 computed 
without crystal field (left), including the crystal field of the covalent type (center), and the crystal field of both covalent and 
Madelung types (right). 



TABLE III: Magnetic interactions (Jrr'), Hartree-Fock energies (Ehf), and total energies (E to t) in the orthorhombic phase of 
YVO3 (T<77 K). The energies are measured from the most stable magnetic state in meV per one formula unit. The magnetic 
interactions are measured in meV. The total energy is defined as Hartree-Fock energy plus correlation energies (Ec) in the 
second order of perturbation theory: E to t=EuF+Ec- 



phase 


J12 


J13 


J24 


J34 


-Ehf 


-Etot 


F 


-1.4 


-3.6 


-3.6 


-1.4 


21.7 


27.1 


A 


-2.5 


-3.7 


-3.7 


-2.5 


14.6 


17.3 


C 


-4.7 


-4.9 


-4.9 


-4.7 


10.1 


11.7 


G 


-4.4 


-4.8 


-4.8 


-4.4 








flip 


-4.0 


-4.8 


-3.9 


-1.8 


11.6 


14.0 



The absolute values of J12 and J13 obtained in the HF calculations for the G-type AFM phase are underestimated 
by about 1 meV in comparison with the experimental parameters extracted from the fit of the spin-wave spectra 
(Ji2=Ji3= — 5.7±0.3 meV). This seems to be reasonable because the HF method is a single-Slater-determinant 
approach, which does not include the correlation effects. The magnitude of the correlations energy depends on the 
magnetic state and is expected to be larger in the case of the AFM spin alignment, where the net magnetization 
is zero and the choice of the many-electron wave function in the form of a single Slater determinant is always an 
approximation^ On the other hand, the saturated ferromagnetic state can be described relatively well by a single 
Slater determinant. All these trends are clearly seen in the total energy calculations, which take into account the 
correlation effect in the second order of perturbation theory (Table UTTjl . The correlations additionally stabilize the 
G-type AFM ground state relative to other magnetic states. However, it does not change the order of the magnetic 
states. 

Unfortunately, it is not possible to estimates the effects of correlations on the interatomic magnetic interactions 
directly, using Eq. (JJJ). Nevertheless, one can try to use the total energy differences between different magnetic states 
and map them onto the Heisenberg models This is a cruder approximation, which implies that the orbital ordering is 
completely quenched by the crystal distortion and does not depend on the magnetic state. We will use it here only in 
order to get a qualitative idea about the impact of correlation effects on the interatomic magnetic interactions. Then, 
the standard HF approximation yields J\2= —3.3 meV and J13— —4.3 meV, the second order perturbation theory for 
the correlation energy yields Ju= —4.1 meV and J13— —5.4 meV. Therefore, the 1 meV difference between parameters 
of magnetic interactions obtained in the HF calculations and the experimental spin-wave dispersion data can be 
naturally attributed to the correlation effects. This example also clearly shows at the importance of more rigorous 
treatment of the correlation effects and necessity to go beyond the single-Slater-determinant HF approximation for 
the tig perovskite oxides. 
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FIG. 12: (Color online) Distribution of the charge density around V-sites in various magnetic phases of monoclinically distorted 
YVO3, as obtained in Hartree-Fock calculations. Different magnetic sublattices are shown by different colors. 



2. High-temperature monoclinic phase (77 K < T < 116 K) 

The monoclinic distortion creates two inequivalent types of V-sites, which lie in different ab-planes: (1,2) and (3,4) 
in Fig. |21 This leads to the new scheme of the i2g-level splitting (Fig. (SJ. Energetically, the new scheme is rather 
similar to the previous one, observed in the orthorhombic phase. In both planes it lowers the energies of two ti g 
levels. The energy gap, which separates the highest level from the middle one is 101 and 128 meV for the sites 1 and 
3, respectively. However, the type of the orbitals which are split off by the distortion is different (Table[IJ|. This leads 
to the new type of stacking between the planes, which reminiscent the G-type orbital ordering. 

However, the orbital ordering is not completely quenched by the crystal distortion and there is a substantial variation 
of the orbital structure depending on the magnetic state, which can be seen already in the distribution of the charge 
density in Fig. 1121 Nevertheless, the basic G-type orbital ordering pattern is clearly seen in all magnetic structures. 
This appears to be sufficient to stabilize the experimentally observed C-type AFM phase, which becomes the ground 
state already in the HF approach (Table llV|l . The correlation effects play a very important role. They additionally 
stabilize the C-type AFM ground state and reverse the order of other magnetic states (e.g., F and A in Table llVl) . 

The orbital ordering in the plane (3,4) clearly reminiscent the one observed in the orthorhombic phase (Fig. I10[l . 
The shape of orbitals in the plane (1,2) appears to be more distorted. 

The behavior of interatomic magnetic interactions in the high-temperature phase of YVO3 has attracted a con- 
siderable attention recently. The experimental spin-wave spectrum shows a clear splitting into acoustic and optical 
branches, which are separated by a 5 meV gap in the middle, q m =(0, 0, 7r/2c), point of the first Brillouin zone along 
the [001] direction^ The splitting has been initially attributed to the dimerization effect associated with an orbital 
Peierls state, which should lead to the alternation of the strong and weak ferromagnetic bonds along the c axis. 16 
However, more recently the puzzling features of the experimental spectra have be naturally explained by the difference 
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TABLE IV: Magnetic interactions (Jrr/), Hartree-Fock energies (Ehf), and total energies (E to t) in the monoclinic phase of 
YV0 3 (77 K < T < 116 K). See Table EH for other notations. 



phase 


J12 


Jl3 


J24 


J34 


-Ehf 


-Etot 


F 


0.1 


2.7 


2.7 


-3.5 


11.7 


17.5 


A 


-0.3 


2.5 


2.5 


-3.4 


14.0 


17.1 


C 


-0.9 


2.2 


2.2 


-4.5 








G 


-1.6 


1.8 


1.8 


-4.8 


6.6 


7.2 


flip-I 


-1.8 


1.6 


2.5 


-3.4 


11.5 


14.2 


fiip-II 


-0.4 


2.3 


2.8 


-4.8 


4.4 


6.5 



of the exchange coupling constants in the adjacent ab-planes, which is expected for the monoclinic C\ h symmetry^ 
This effect is clearly seen in our HF calculations: while the AFM exchange coupling in the plane (3,4) does not change 
so much in comparison with the orthorhombic phase, the one in the plane (1,2) drops by almost 4 meV (referring to 
the C-type AFM state in Table Hy), The value of the gap in the point q m can be estimated in the linear spin-wave 
theory as Agw=2Ji3| yl — 4J12/J13— y'l — 4J 3 4/Ji 3 |. Then, using results of HF calculations we obtain Asw~6-2 
meV, which is in fair agreement with the experimental finding. We can further speculate that the ferromagnetic 
coupling J 13 is overestimated in the HF approximation due to the lack of the correlations effects (in the next section 
we shall sec that this is indeed the case for the totally ferromagnetic YTi0 3 ). Therefore, the correlation effects will 
probably yield a smaller value of AgW: which is proportional to Ji 3 . 



B. YTiOa 

YTi0 3 is a d 1 compound. The lattice distortion splits off one t 2g -\eve\ to the low-energy part of the spectrum 
(Fig. ^J. This is just enough for trapping one d electron at each Ti-site. Therefore, the situation is similar to YV0 3 . 
The lowest ^g-levcl is separated from the middle one by a 109 meV gap, which is comparable with the magnitude of 
the CF splitting in YV0 3 . 

The lattice distortion stabilizes the orbital ordering, which is shown in Fig. In this case, the orbital ordering 
is strongly quenched by the distortion and not only the charge density but also the parameters of the magnetic 
interactions (TableEJ, which are sensitive to the orbital ordering, only weakly depend on the magnetic state. Partly, 



TABLE V: Magnetic interactions (Jrr'), Hartree-Fock energies (Ehf), total energies (E to t), and superexchange energies (Ese) 
in YTiOa. Ese is defined as the total energy in the superexchange approximation. Other definitions arc explained in Tabic UTTI 



phase 


J12 


Jl3 


J24 


J34 


Ehf 


Etot 


Ese 


F 


3.9 


1.2 


1.2 


3.9 











A 


3.9 


1.1 


1.1 


3.9 


2.0 


0.8 


0.1 


C 


3.2 


1.1 


1.1 


3.2 


14.4 


10.9 


12.1 


G 


3.2 


1.0 


1.0 


3.2 


16.2 


12.5 


13.9 


flip 


3.2 


1.1 


1.1 


3.9 


8.2 


6.1 


6.7 



this is because YTi0 3 has the largest on-site Coulomb interaction U amongst considered perovskite oxides ( Tabic ITI) . 
Therefore, the superexchange contribution to the orbital ordering is expected to be smaller in comparison with the CF 
splitting. The obtained orbital ordering is in the excellent agreement with results of LDA+£7 calculations by Sawada 
and Terakura^i and the experimental measurement using the nuclear magnetic resonance (Ref. |jj), the polarized 
neutron diffraction (Ref. 8), and the soft x-ray linear dichroism (Ref. ITcJl . 

The observed orbital ordering patter is compatible with the ferromagnetic ground state, which easily emerges at the 
level of HF theories ,iiL2Zi The same trend is clearly seen in our calculations, and both the order F^A^C^flip^G and 
the total-energies difference between different magnetic states obtained in the HF approximation are well consistent 
with the LDA+C7 results of Sawada and Terakura^ 

The correlation effects are important. Similar to YV0 3 , they tend to additionally stabilize the AFM configurations 
and destabilize the ferromagnetic ground state. The situation is very fragile, and after taking into account the 
correlation effects, the energy difference between F state and the next A- type AFM state becomes only 0.1-0.8 meV 
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FIG. 13: (Color online) Distribution of the charge density around Ti-sites in various magnetic phases of YTiOa, as obtained 
in Hartree-Fock calculations. Different magnetic sublattices are shown by different colors. 



per one formula unit. For the d 1 compounds, we can apply two independent schemes for studying the correlation 
effects: one is the second order perturbation theory and the other one is the theory of superexchange interactions, 
which takes into account the multiplet structure of the excited states. For YTi03, both scheme yields very consistent 
results, apart from the small quantitative difference, which is inevitable for two different approximations. 

Yet, the magnetic behavior of YTiC>3 poses several open questions, which are not fully understood. 

First of all, YTiC>3 is an exceptional example amongst ti g perovskite oxides, because the ferromagnetic ground state 
can be anticipated already on the basis of the canonical Goodenough-Kanamori- Anderson rules for the superexchange 
interactions in the simple cubic structure. This immediately revives the idea of Kugel and Khomskii about the 
superexchange-driven orbital ordering and the concomitant Jahn- Teller distortion, which was intensively discussed in 
the context of KCuF3^i Then, one may ask whether the experimental orbital ordering in YTi03 can be stabilized by 
the pure superexchange mechanism, without appealing to the CF splitting. This can be easily checked by substituting 
/i^=0 into kinetic-energy part of the model Hamiltonian. Surprisingly, the orbital ordering in the ferromagnetic state 
is practically the same with and without the CF splitting fRef. I14fl . The interatomic magnetic interactions, Ji2=4.3 
meV and Ji3=1.6 meV, are also consistent with the data listed in Table [V] and which include the effects of the CF 
splitting. This naturally explain results of our previous work fRef l28l) . where similar orbital ordering and interatomic 
magnetic interactions have been obtained without the Madelung term in the CF splitting. Therefore, it is tempting 
to conclude that the orbital ordering in YTi03 is driven by the superexchange interactions, and the lattice distortion 
simply follows the anisotropic distribution of the charge density associated with this orbital state. However, the 
situation is not so simple, because our calculations have been performed in the room temperature structure (T=293 
K, Ref. 0), which is much higher than the Curie temperature (Tc«30 K, Ref. 0). Furthermore, the orbital ordering 
shown in Fig. ^] in the absence of the CF splitting takes places only in the ferromagnetic phase: had we changed the 
magnetic state, our orbital ordering would have been also different. Therefore, a more plausible scenario for YTi03 is 




FIG. 14: (Color online) Orbital ordering in the ferromagnetic phase of YTi03 computed with (left) and without (right) 
crystal-field splitting. 



that the lattice distortion goes first and sets up a particular form of the CF splitting, which stabilizes the experimental 
orbital ordering pattern and the ferromagnetic ground state. Nevertheless, the good agreement between two orbital 
states shown in Fig. 1141 is really curious. Is it a simple coincidence or there is some physical meaning behind this 
result? We would like to emphasize again that the situation is totally different from YVO3 shown in Fig. 1111 

Another group of questions is related with the behavior of interatomic magnetic interactions. The first puzzling 
feature is the nearly isotropic experimental spin-wave spectrum reported by Ulrich et al^ which cannot be explained 
in terms interatomic magnetic interactions derived from the first-principles electronic structure calculations. For 
example, the exchange parameters listed in Table [V] are clearly anisotropic, and the interactions along the c-axis are 
much weaker than in the ab-plane. A similar conclusion is expected from the analysis of the total-energy differences 
reported by Sawada and Terakura?22, Since the anisotropy of magnetic interactions is directly related with the particular 
form of the orbital ordering, it seems that the inelastic neutron-scattering data by Ulrich et al. are in an apparent 
disagreement not only with the results of the first-principles electronic structure calculations but also with a number of 
other experimental data, which report the same type of the orbital orderingA&ifi Therefore, what is so special about 
the inelastic neutron-scattering measurements and why different experimental methods probing the orbital state yield 
qualitatively different conclusions in the case of YTiC>3? 

The behavior of interatomic magnetic interactions predetermines not only in the form of the spin-wave spectrum, 
but also in the absolute value of To If the magnetic properties of YTiC>3 are indeed controlled by the large lattice 
distortion, which is set up far above To it should to be a good Heisenberg ferromagnet. This is directly seen in our 
HF calculations, where the parameters of interatomic magnetic interactions practically do not depend on the magnetic 
state (Table lVjl . Then, the applicability of the Heisenberg model is no longer restricted by infinitesimal rotations of the 
spin magnetic moments, and Tc can be easily evaluated using the standard expressions, which are well known from the 
theory of Heisenberg magnets^ The simplest one is the mean-field formula: 3fceT c F —A J12+2 J13, where the prefactor 
S(S+1) is already included in the definition of our exchange parameters, although with some approximations for the 
spin l/2i«£ By combining this expression with the HF approximation for the exchange interactions, one finds T C F =64 
K, which exceeds the experimental value by factor two. However, T c F does not include spontaneous fluctuations and 
correlations between the motion of the neighboring spins. This is exactly the point where the anisotropy of exchange 
interactions can help to reduce the theoretical value of Tc- Indeed, according to the Mermin- Wagner theorem, 58 
the two-dimensional Heisenberg model does not support any long-range spin order at any nonzero temperature. 
Therefore, since for Ji3/Ji2 < Cl the system will eventually approach the two-dimensional limit, it is reasonable to 
expect a substantial reduction of Tc- In order to describe these effects quantitatively, one can use the spherical 
approximation for the Heisenberg model^ according to which 3fcBTc=l/Ek(^o~4) _1 ], Jk=J2n' ^RR' e ' k R j an d 
the summation over k is restricted by the first Brillouin zone. Then, using parameters extracted from the HF 
calculations we obtain Tc=36 K, which can be further reduced by taking into account the correlation effects. For 
example, by mapping the total energies obtained in the second order of perturbation theory and in the theory of 
superexchange interactions onto the Heisenberg model, we obtain Tc= 27 and 28 K, respectively, which is in fair 
agreement with the experimental data. 
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FIG. 15: (Color online) Distribution of the charge density around V-sites in various magnetic phases of LaVOs, as obtained 
in Hartree-Fock calculations . Different magnetic sublattices are shown by different colors. 



C. LaV0 3 

The monoclinic LaVC>3 has two inequivalent V-sites. Contrary to YVO3, these sites differ not only by the direction, 
but also by the magnitude of the CF splitting (Table 0, which is 78 and 152 meV for the sites 1 and 3, respectively 
(referring to the splitting between middle and highest t2 g levels). Therefore, already from this very simple analysis 
of the CF splitting one can expect very different behavior of the orbital degrees of freedom in different ab-planes: 
the strong quenching in the plane (3,4), and a relative flexibility in the plane (1,2) (Fig. [2J. Thus, the situation is 
qualitatively different from YVO3. 

These trends are clearly seen in the HF calculations (Fig. ll5|) : the orbital ordering in the plane (1,2) strongly depend 
on the magnetic state and one can clearly distinguish two types of the orbital-ordering pattern realized, on the one 
hand, in the states F and C, and, on the other hand, in the states A and G. Therefore, it is perhaps right to say 
that in LaV03, the experimental orbital ordering is partly driven by the magnetic ordering via the Kugel-Khomskii 
mechanism, - 1 This seems to agree with the experimental data, which show that in La-based compounds, the orbital 
ordering develops few degrees below the magnetic Neel temperature (Tn)iS This is again different from YVO3, for 
which the orbital-ordering temperature is substantially higher than Tn- 

The change of the orbital ordering is reflected in the behavior of interatomic magnetic interactions, which not only 
depend on the magnetic states, but can even change the signs (Table IVlj) . In such a situation, the total energy may 
have several local minima, realized for those magnetic states where the signs of interatomic magnetic interactions are 
consistent with the type of the imposed spin ordering. We have found at least two such minima, corresponding to the 
C-type AFM ground state and the G-type AFM state, which has higher energy. Thus, we do not quite agree with the 
conclusions about the complete quenching of the orbital ordering in these distorted perovskite compounds. 29 This is 
not generally true and LaV03 is clearly an exception. However, the CF splitting will also prevent the formation of the 
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TABLE VI: Magnetic interactions (Jrr/), Hartree-Fock energies (Ehf), and total energies (E to t) in LaVOs. See Table HTT1 for 
other notations. 



phase 


Jl2 


Jl3 


J24 


J34 


-Ehf 


-E to t 


F 


-5.1 


6.6 


6.6 


-1.7 


21.0 


30.8 


A 


3.8 


2.1 


-4.5 


-2.4 


20.6 


24.1 


C 


-4.8 


6.0 


6.0 


-1.8 








G 


-6.3 


-4.4 


-4.4 


-2.4 


7.6 


11.1 


flip-I 


7.7 


0.6 


5.8 


-1.7 


9.8 


15.0 


flip-II 


-6.2 


-4.0 


5.8 


-3.1 


12.7 


16.9 



orbital singlet states, which were employed in order to explain the appearance of the C-type AFM antiferromagnetism 
in the theory of orbital fluctuations iSi Actually, such a singlet state conflicts with the C\ h symmetry of the monoclinic 
phase. 

Since the interatomic magnetic interactions depend on the magnetic state, the simple Heisenberg model may be 
used only for the analysis of local perturbation around each magnetic state. Then, it is reasonable to expect a gap 
Asw~6.8 meV between acoustic and optical branches of the spin-wave spectrum, similar to the one observed in the 
C-type AFM phase of YVO31S& However, the Heisenberg model is no longer valid for the analysis of the transition 
temperature (unlike in YTiOa), which should take into account a possible change of the the orbital states in the 
course of thermodynamic average. 

Similar to YVO3 and YTiC>3, the correlation effects play a very important role also in LaVC>3 and additionally 
stabilize the C-type AFM ground state. 

D. LaTi0 3 

LaTiOa is a puzzling system. It has the smallest CF splitting (Fig. [SJ among the distorted perovskite oxides, which 
formally leaves a room for the orbital fluctuations. On the other hand, the possible variation of the orbital order 
appears to be bounded by certain constraint conditions. For example, although the orbital ordering depends on the 
magnetic state, this dependence is not particularly large, as it is clearly seen from the HF calculations, where the basic 
shape of the orbital-ordering pattern remains the same for different magnetic states (Fig. I16H . There are certainly 
some variations of the orbital ordering, which can be seen already on the plot. However, these variations do not seem 
to change the qualitative conclusion about the form of interatomic magnetic interactions and the magnetic ground 
state of LaTi03. 

Unfortunately, this conclusion is not consistent with the experimental data. In this sense, there is a clear difference 
of LaTi03 from other perovskite compounds considered in this work, despite the fact that we have used absolutely 
the same procedure for construction and solution of the model Hamiltonian. 

The magnetic ground state is expected to be of the A-type, as it is clearly seen from the total-energy calculations as 
well as from the behavior of interatomic magnetic interactions ( Table IVlIj l. although experimentally, it is totally an- 
tifcrromagnctic G-typeM> The magnetic interactions are sensitive to the magnetic ordering. However, other magnetic 
states appear to be unstable and the form of interatomic magnetic interactions in each magnetic state systematically 
leads to the A-type antiferromagnetism. This conclusion is totally consistent with results of our previous workj2& 
which neglects the Madelung contribution to the CF splitting. 



TABLE VII: Magnetic interactions (Jim/), Hartree-Fock energies (-Ehf), total energies (E to t), and superexchange energies 
(-Ese) in LaTiC>3. See Tables HTT1 and IV1 for other notations. 



phase 


J12 


J13 


J24 


J34 


Ehf 


Etot 


E SE 


F 


4.5 


-1.2 


-1.2 


4.5 


5.0 


17.7 


1.6 


A 


3.6 


-3.3 


-3.3 


3.6 











C 


1.0 


-2.0 


-2.0 


3.4 


19.6 


26.3 


14.2 


G 


2.0 


-4.9 


-4.9 


2.0 


11.5 


11.0 


9.0 


flip 


1.3 


-4.5 


-0.4 


4.6 


7.7 


11.4 


5.7 



Then, what is missing in our calculations - or maybe even more generally - in our understanding of the magnetic 
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FIG. 16: (Color online) Distribution of the charge density around Ti-sites in various magnetic phases of LaTiOs, as obtained 
in Hartree-Fock calculations. Different magnetic sublattices are shown by different colors. 



properties of LaTiC>3? Below we discuss several plausible scenarios. 

(i) One possibility is that the effect of the crystal distortion maybe still underestimated. Particularly, we tried to 
follow the idea of Refs.|2^ and0, an d additionally scaled the contribution of the Madelung term to the CF splitting 
by multiplying the right-hand side of Eq. {2} by a constant. This corresponds to the change of the dielectric constant, 
which was treated as an adjustable parameter in Refs. |2^ and |3 We have found that in order to obtain the 
experimental G-type antiferromagnetic ground state, the dielectric constant should be reduced by a factor four or 
five (Fig. IT?}) . Then, the exchange interactions become nearly isotropic (Ji2=— 3.3 and J13—— 3.5 meV) and do not 
depend on the magnetic state. Hence, LaTiOs is expected to be a good Heisenberg antiferromagnet, in agreement 
with the experimental inelastic neutron-scattering data. The latter reveals nearly isotropic behavior of the spin- wave 
spectrum with J12 = J13 = — 15. 5S" 2 « —3.9 meVM Thus, we are totally agree with the authors of Refs. |2^ and 
|2J that the Madelung term alone could explain many experimental features of LaTiC>3. The only problem is that, 
according to the electronic structure calculations, the effect is too small. One can of course try to blame LDA for 
this failure. However, why does this problem occur only for LaTiOs while for other compounds our method works 
reasonably well? We believe that if the story about the structural origin of the G-type antiferromagnetism in LaTi03 
make a sense, it is more likely that the real magnitude of the crystal distortion in LaTiOs may be still undisclosed 
experimentally This seems to be reasonable, because the structural data for the distorted perovskite oxides are still 
in the process of steady refinement .iSiiiiSi 

(ii) Other scenarios are related with the correlation effects, which are not included in the HF calculations. There is 
no doubts that they must play an important role in LaTiOs. However, it seems that there is no simple scheme which 
would allow to include these effects in the electronic structure calculations. The second-order perturbation theory 
and the theory of superexchange interactions, which we have tried, are definitely not enough. They do substantially 
change the total energies of the HF method. However, the conclusions strongly depend on the approximation which 
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FIG. 17: (Color online) Total energies measured relative to the G-type antilerromagnetic state obtained in Hartree-Fock 
calculations with exaggerated crystal-field splitting of the Madelung type. The latter as been multiplied by the scaling factor 
shown in the absciss axis. 

we use (Table fVIIfl . For example, from the second order perturbation theory, it seems to be clear that the correlation 
effects tend to stabilize the G-type AFM state: the correlations change the order of the magnetic states and somewhat 
lower the energy of the G-type AFM state relative to the A-state. The latter trend is also seen the superexchange 
approach. However, the superexchange method tends to lower also the energies of other magnetic states, apparently 
through the small change of the orbital ordering. This is in straight contrast with the more distorted YTi03, were 
two different methods provided rather consistent explanation for the role played by the correlation effects (Table Ivjl . 
(iii) The theory of orbital liquid is just an opposite case to the theory of CF splitting as these two effects are 
incompatible with each other. Although the formation of the orbital liquid in the cubic ti g lattice is a many-electron 
effect, the necessary prerequisite, which should exist already at the mean-field level, is an infinite degeneracy of the 
magnetic ground stated Although our CF splitting for LaTi03 is small, that is sometimes regarded as the strong 
support for the orbital liquid theory, we do not observe such a degeneracy of the Hartree-Fock ground state: all 
HF calculations steadily converge to a single solution for the orbital-ordering pattern, which is shown in Fig. ^1 
irrespectively on the starting conditions and the size of the supercell. Therefore, although the correlation effects 
beyond the HF approximation are certainly important in the case of LaTiOs, it does not necessary mean that the 
ultimate result of electron correlations must be formation of the orbital liquid state. We hope that even in LaTi03, 
the correlation effects can be systematically included using the regular perturbation theory expansion around the 
nondegenerate HF ground state, while the second order of this expansion is simply not enough. In such a situation, 
it is perhaps more practical to derive an approximate ground state from the diagonalization of a many-electron 
Hamiltonian matrix constructed in the basis of some limited number of specially selected Slater determinants, as it is 
done for example in the path-integral renormalization group method^S 

Thus, one of the challenging problems in the theory of distorted perovskite oxides remains to be the explanation 
of the G-type AFM ground state in LaTiOs. Definitely, we need a more rigorous theory for the correlation effects. 
But, will it be enough, or do we need a more radical refinement of our starting model, given by Eq. QJ, in the case 
of LaTi0 3 ? 



E. Spin-Orbit Interaction and Magnetic Ground State 

The spin-orbit interaction in the distorted perovskite oxides will generally leads to the noncollinear spin arangement, 
which obeys certain symmetry rules*' 1 The spin magnetic moments, aligned along one of the orthorhombic axes, will 
be subjected to certain rotational forces, coming both from the Dzyaloshinsky-Moriya interactions and from the 
single-ion anisotropy energies^S which lead to the spin canting and the appearance of nonvanishing components of 
the spin-magnetization density along two other directions. The type of the magnetic ordering for all three projection 
of the spin-magnetization density is generally different. Thus, each magnetic structure can be generally abbreviated 
as X-Y-Z, where X, Y, and Z is the type of the magnetic ordering (F, A, C, or G) formed by the projections of 
the spin magnetic moments onto the orthorhombic axes a, b, and c, respectively. The orbital magnetic structure has 
the same symmetry, although it has a different origin of the canting, which comes mainly from the interplay of the 
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spin-orbit interaction with the CF splitting at each transition site. The spin and orbital magnetic moments are not 
generally collinear to each otheriSi 

Results of HF calculations, which take into account the spin-orbit interaction, are summarized in Table IVIIII 
The magnetic ground state of YTiC>3 is G-A-F, in agreement with the neutron scattering data. 9 The ferromagnetic 

TABLE VIII: The type of the magnetic ground state, the vectors of the spin ((is) and orbital (pl) magnetic moments (in /xs), 
and the values of the band gap (E g , in eV) obtained in Hartree-Fock calculations. For the magnetic ground state, three capital 
letters denote the type of the magnetic ordering for three projections of the magnetic moments onto the orthorhombic axes 
a, b, and c, respectively. For the spin and orbital magnetic moments at the transition-metal sites, these three projections are 
specified by the vectors (is and fj,L, respectively. The positions of the transition-metal sites are shown in Fig. [5] 



compound 


phase 


site 


ground state 






Ms 








Ml 




E g 


YT1O3 





1 


G-A-F 


(- 


-0.00,- 


-0.09, 


0.84) 


(- 


-0.05, 


0.01,- 


-0.03) 


1.1 


LaTi0 3 





1 


GF-A 


( 


0.02,- 


-0.17, 


0.76) 


( 


0.10, 


0.03,- 


-0.08) 


0.6 


YVO3 





1 


F-GG 


(- 


-0.02, 


0.00, 


1.65) 


( 


0.00,- 


-0.00,- 


-0.17) 


1.2 


YVO3 


m 


1 


GA-C 


(- 


-0.74, 


0.08, 


1.48) 


( 


0.07,- 


-0.04,- 


-0.16) 


1.0 






3 




(- 


-0.78,- 


-0.03, 


1.48) 


( 


0.05, 


0.04,- 


-0.07) 




LaV0 3 
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moment is aligned along the c-axis, and the canting angle is relatively small. The absolute values of the spin and 
orbital magnetic moments are different from those reported in our previous work^S The difference is related with the 
additional transformation to the LMTO basis, which allows to calculate the local magnetic moments directly at the 
transition-metal sites, whereas in Ref. |2S| they have been calculated in the Wannier basis, which had a substantial 
weight at the oxygen sites. This comparison clearly shows that the covalent effects play a very important role and 
allow to explain a substantial reduction of the local magnetic moments at the transition-metal sites. As it was already 
discussed in Sec. IV Bl the correlation effects in YTiC>3 favor the AFM coupling and systematically lower the energies 
of all AFM states relative to the ferromagnetic one (see Table 0. After including the spin-orbit interaction, this will 
lead to the additional spin canting away from the collinear ferromagnetic arrangement. For example, in the framework 
of superexchange approach we have obtained the following values of the spin and orbital magnetic moments (in hb, 
referred to the site 1): ns=(— 0-02, —0.29, 0.77) and fiL=(— 0.05, 0.02, —0.03). It readily explains the experimental 
values reported for the F and G components of the magnetic moments, correspondingly along the c and a-directions^ 
while the A-type AFM component along the b-direction is somewhat larger in our calculations. The reason is not 
completely clear. The orbital magnetic moments in YTiOs are strongly quenched by the crystal field. 

The orthorhombic YVO3 has nearly-collinear magnetic structure, where the G-type AFM moment is aligned along 
the c-axis. The total magnetic moment (/is+Mi)ll c= l-48/is parallel to the c-axis is in excellent agreement with the 
experimental value of 1.45/zb reported by Blake et al. and corresponding to T=65 K. 14 Ulrich et al. reported somewhat 
larger value of 1.72yUs, also oriented along c. The weak ferromagnetic component along the a-direction has been also 
observed experimentally. 13 Details of magnetic ordering in the monoclinic phase of YVO3 are somewhat controversial. 
Blake et al. reported the C-type AFM ordering for both b- and c-components of the magnetic moments in the 
orthorhombic Pbnm notations, which correspond to the a- and c-directions in the monoclinic P1\ja notations^ 
This is totally consistent with our finding. The quantitative difference can be naturally explained by the finite 
temperature effects in the intermediate phased Furthermore, we predict the A-type antiferrimagnetic ordering for 
the remaining b-component (in the P2\/a notations), which implies nearly antiferromagnetic alignment of the b- 
projections of the magnetic moments in the adjacent ab-planes. However, since the sites 1 and 3 are not fully 
equivalent in the monoclinic structure, the b-projections do not compensate each other, and the system exposes the 
net magnetic moment, which can couple to the magnetic fields The antiferrimagnetic ordering along the b-axis 
is also consistent with the temperature-induced magnetization reversal behavior of YVO3 observed by Ren^ and 
could be a natural explanation for this effect. More generally, it is right to say that the magnetic coupling along the 
c-dircction of the monoclinic phase is always cither ferromagnetic or antiferrimagnetic, because the sites 1 and 3 are 
not equivalent. 

The C-A-C magnetic ground state appears to be different from the magnetic structure reported by Ulrich et al., 
who have observed the C-type AFM ordering for the a- and b-components, and the G-type AFM ordering for the 
remaining c-componcnt (apparently, in the Pbnm notations)^ One possible explanation for this difference could be 
the coexistence of several magnetic states in a narrow energy range. Such a behavior has been indeed observed in 
our HF calculations: in addition to the C-A-C state we could obtain another self-consistent solution corresponding 
to the A-C-F phase with a slightly higher energy (about 0.04 meV per one formula unit, which is of the order of the 
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magnetocrystalline anisotropy energy). The new phase has the following magnetic moments (in fiB, where the first 
and second lines corresponds to the sites 1 and 3, respectively): 

/4 = (-0.07, 1.65, 0.05), A*i = ( 0.05,-0.15,-0.02), 
A*! = ( 0.07, 1.67, 0.03), n 3 L = (-0.09,-0.07,-0.02). 

Although the exact form of the magnetic state A-C-F is still different from the observation by Ulrich et al.^ we can 
speculate that they probably used a different experimental setup which yielded the realization of another magnetic 
phase, which was different from the finding by Blake et al. fRef. Ibfl) and Ren et al. fRef. Il3j) . At least, the C-type 
AFM component parallel to the ab-plane is qualitatively consistent with the form of the A-C-F phase obtain in our 
HF conclusions, and to certain extent the results of HF calculations can be further modified by the correlation effects. 

In the case of LaVC>3 we could find three stable solutions: A-C-A, C-A-C, and C-F-G. The first two are similar 
to the A-C-F and C-A-C states, emerging in monoclinic YVO3. The only difference is that the A-C-A state has 
lower energy and, therefore, is realized as the ground state in the case of LaV03. Like in YVO3, the A-type state 
corresponds to the antiferrimagnetic ordering. Therefore, LaVC>3 is expected to have the net magnetic moment in 
the ac-plane. The third (C-F-G) solution has nearly collinear spin structure: 

M L = (-0.01,-0.06, 1.60), nl = { °- 01 > 0.13,-0.34), 
M | = (-0.03, -0.02, -1.61), nl = ( 0.01, 0.04, 0.14). 

It corresponds to the stable G-type AFM phase emerging without the spin-orbit interaction (see Table ITTTf) . Never- 
theless, the A-C-A state appears to be well separated from the C-A-C and C-F-G states, correspondingly by 0.19 and 
6.50 meV per one formula unit. 

The quenching of the orbital magnetic moments in two different sublattices of monoclinic YVO3 and LaV03 nicely 
correlates with the values of the CF splitting (Fig.|SJ|, where larger CF splitting at the site 3 results in smaller orbital 
magnetization. 

It is probably meaningless to discuss the relativistic effects in LaTi03, where we could not reproduce the correct 
magnetic ground state. We could agree with the criticism risen by Haverkort et al. (Ref. l4lTl that our CF splitting 
alone does not explain details of their spin-resolved photoemission spectra (actually, our value of the parameter (L • S) , 
obtained in the HF calculations after the transformation to the LMTO basis and radial integration over the Ti-sphcre 
is about —0.13, which exceeds the experimental value by factor two). However, it does not make a sense to present 
as an alternative the results of calculations yielding the same A-type AFM ground stated which totally agrees with 
our finding and (unfortunately) disagrees with the experiment, even though these calculations yield somewhat larger 
value of the CF splitting for LaTiOs. As it has been already discussed in Sec. lIII Al different values of the CF splitting 
obtained by different authors are most probably related with the nonunique choice of the Wannier functions for the 
t2 g bands of the distorted perovskite oxides. One can formally adjust the theoretical CF splitting in order to meet 
certain demands of some particular class of the experimental data. However, will it solve a more fundamental problem 
of the magnetic ground state of LaTiOs? 

In addition, we show in Table IVlTll the values of the band gap obtained in the HF calculations for the magnetic 
ground state. For YTi03, YVO3, and LaV03 there is a good agreement with the experimental optical dataii**^ 
However, for LaTiOs the experimental gap is much smaller (~0.1 eV)& This may indicate again at the particular 
importance of correlation effects in LaTiOs. 

VI. SUMMARY AND CONCLUDING REMARKS 

The main purpose of this work was to make a bridge between first-principles electronic structure calculations 
and model approaches for the strongly-correlated systems, and illustrate how it works for the series of distorted 
t2g perovskite oxides. The whole plan included three major steps: first-principles electronic structure calculations 
— > construction of the model Hamiltonian for the isolated ti g bands near the Fermi level — > solution of this model 
Hamiltonian using several different approaches. The choice of the distorted t2 g perovskite oxides was motivated by 
the fact that they represent a good example of systems where it is practically impossible to construct a relevant 
model without the impact from the first-principles calculations: simply, the lattice distortion is too complex and 
there are too many model parameters, which cannot be fixed in unbiased way. In this sense, we strongly believe that 
any theoretical model for such complex oxide materials should be based on the results of first-principles electronic 
structure calculations. Otherwise it could be just an abstract mathematical construction deprived of clear physical 
grounds. 

The present work clearly demonstrates that nowadays the idea of constructing the parameter-free model Hamiltoni- 
ans for the strongly-correlated systems is quite feasible: all model parameters in our work, including the intraatomic 



28 



Coulomb interactions, have been derived from the first-principles calculations using the method proposed in Ref. l34l 
Apart from the approximations inherent to the method, the procedure of constructing the model Hamiltonian was 
totally parameter- free, and our analysis did not rely on the use of any adjustable parameters. Therefore, it is remark- 
able that using such a parameter-free approach we could propose a consistent explanation for a number of puzzling 
properties of the distorted t 2g perovskite oxides. These first results are really encouraging and we would like to 
hope that the method proposed in Ref. 03 can be successfully applied in the future for the analysis of electronic and 
magnetic properties of other narrow-band compounds. ^^^^ 

At the mean-field Hartree-Fock level, the results of all-electron full-potential LDA+J7 calculations (Refs. I27l29() 
can be successfully reproduced in our model approach for the isolated t 2g bands. We could easily rationalize the 
main results of these relatively heavy calculations and elucidate the main microscopic interactions responsible for the 
formation of different magnetic structures in the case of YTiC>3, YVO3, and LaV03. We argue that the nonsphericity 
of the Madelung potential should be an indispensable ingredient of electronic structure calculations, and the results 
of commonly used atomic-spheres-approximation should be corrected in order to include these effects. 

The crystal distortion plays an important role in the physics of t 2g perovskite oxides. At least for YTiC>3, YVO3, 
and LaVOs , the knowledge of the experimental lattice parameters and the atomic positions greatly helps in explaining 
the magnetic properties of these compounds. Of course, some questions still remain. Particularly, what is the origin of 
this distortion? Why is it so different for different compounds, that is finally manifested in the formation of completely 
different magnetic structures? Very similar arguments have been used for LaMnOs, which is a characteristic example 
of the Jahn- Teller distorted e g perovskite oxides. Particularly, it was argued that the experimental distortion not 
only stabilize the A- type AFM ground state of LaMnOsfS&ffiiSLS but is also responsible for the opening of the band 
gapiSS*S£ However, the direction of the Jahn- Teller distortion in LaMnC>3 can be naturally understood in terms of 
anharmonicity of the electron-lattice interactioni^iiSLZii In this sense, LaMnOs is an easy example. Then, is it possible 
to rationalize the behavior of t 2g perovskite oxides in a similar way and come up with some suitable electron-lattice 
model, which would explain not only the direction of the lattice distortion in each particular compound, but also the 
difference between these compounds? 

Finally, we emphasize the importance of correlation effects in the t 2g band of distorted perovskite oxides. Although 
the mean-field Hartree-Fock approach provides a satisfactory description for the magnetic properties of YTi03, YVO3, 
and LaV03, the inclusion of the correlation effects systematically improves the agreement with the experimental data 
for all three compounds. Definitely, LaTiOs is an exceptional example for which we could not obtain the correct G-type 
AFM ground state neither at the level of Hartree-Fock approach nor after including some correlation effects, though 
in a very approximate form. However, we expect that the situation may still change by systematically employing 
more rigorous theories for the correlation effects. The approximations considered in the present work were simply not 
enough in the case of LaTi03. 
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